Excess Kondo resonance in a quantum dot device with normal and superconducting 
leads: the physics of Andreev-normal co-tunneling 



Qing-feng Sun 1 , Hong Guo 1 , and Tsung-han Lin 2 
1 Center for the Physics of Materials and Department of Physics, McGill University, Montreal, PQ, Canada H3A 2T8. 
2 State Key Laboratory for Mesoscopic Physics and Department of Physics, Peking University, Beijing 100871, China 



We report on a novel Kondo phenomenon of interacting 
quantum dots coupled asymmetrically to a normal and a su- 
perconducting lead. The effects of intradot Coulomb interac- 
tion and Andreev tunneling give rise to Andreev bound reso- 
nances. As a result, a new type of co-tunneling process which 
we term Andreev-normal co-tunneling, is predicted. At low 
temperatures, coherent superposition of these co-tunneling 
processes induces a Kondo effect in which Cooper pairs di- 
rectly participate formation of a spin singlet, leading to four 
Kondo resonance peaks in the local density of states, and en- 
hancing the tunneling current. 
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The Kondo effect is a prototypical many-body corre- 
lation effect involving interactions between a localized 
spin and free electrons.! Its regent observation in semi- 
conductor quantum dots (QD)0 □ has generated a great 
deal of theoretical and experimental interests and pro- 
vided rich understanding to many-body phenomena at 
the mesoscopic scaleJj For a QD coupled to two normal 
(N) leads, the-physical origin of the Kondo effect is now 
understoodtoj. Consider a single spin degenerate level 
e c i of the QD such that 6d < [J-n < e d + U , where /xjy is 
the chemical potential of the leads and U the on-site e-e 
interaction energy. An electron of either spin up or spin 
down which occupies td cannot tunnel out of the QD 
because td < AW as well, an electron outside the QD 
cannot tunnel into it unless the on-site Coulomb energy 
U is overcome. Therefore, the first-order tunneling pro- 
cess is Coulomb blockaded. However, due to Heisenberg 
uncertainty, the victuaJ higher-order co-tunneling events 
can still take placairclu by which the electron inside QD 
tunnels out followed by an electron with opposite spin 
tunneling into the QD, on a time scale ~ Ti/\hn — | - 
As a consequence, the local spin is flipped. At low tem- 
peratures, the coherent superposition of all possible co- 
tunncling events gives rise to the Kondo effect in which 
the time-averaged spin in the QD is zero due to frequent 
spin flips: the whole system, QD plus leads, forms a spin 
singlet, and a very narrow Kondo peak located at //jv 
arises in the local density of states (LDOS). 

When one of the leads is a superconductor (S), an- 
other transport process — Andreev tunneling, will occur 
in the normal-superconductor interface in which an in- 
cident electron from the normal side is reflected as a 
hole while a Cooper pair is created in the supercon- 
ductor. Andreev tunneling is very important!] because 



it determines transport properties of many mesoscopic 
superconducting-normal hybrid devices. It is therefore 
not surprising that the Kondo effect in N-QD-S hybrid 
systems has attracted considerable attentionQlla. So far, 
the focus on Kondo effects in N-QD-S devices has been on 
enhancement or reductiasifOf-conductance as compared to 
that of N-QD-N systems0fli!3; and the emergence of sub- 
Kondo peaks in LDOS at ±ALra where A is the super- 
conductor gap energy. In these studies Andreev tunnel- 
ing precesses, in essence, happen alone while the Cooper 
pairs do not participate the formation of spin-singlet. 
However, since the Kondo effect in a QD results from 
co-tunneling processes, for a N-QD-S hybrid system, it is 
very natural to ask: are there co-tunneling processes con- 
sisting of one virtual Andreev tunneling and one virtual 
normal electron tunneling? If there are, can coherent su- 
perpositions of these Andreev-normal co-tunneling give 
rise to a Kondo effect? What are the consequences and 
characteristics of the Kondo effect induced this way? 

It is the purpose of this letter to report our theoreti- 
cal investigation on these issues. In contrast to previous 
worUlniil, we emphasize the possibility of virtual Andreev 
tunneling directly participating the co-tunneling process 
so that the Cooper pairs directly participate formation 
of the spin-singlet: these physical processes give rise to 
Kondo effect in the first place. Our results predicts a 
new co-tunneling process formed by an Andreev tunnel- 
ing and a normal tunneling, and the superposition of this 
type processes induce four Kondo peaks in the LOOS. 

We consider the standard model Hamiltoniant 2 ] of a 
N-QD-S system H = H N + H s + H QD + H T where 
h n = J2ka e Nka{ r7 a klT and H s = J2k* e Skb{ a b ka + 
J2k (Afrfci&-fcT + H.c.) describe the normal lead and the 
superconducting lead, respectively, in here we have set 
[is = 0. Hqd — ^u^dd^dc + Udtd^d]di models the 
QD with a single level having spin index a =f , j and in- 
tradot e-e Coulomb interaction U; Ht = J2kA^ r Na ka d a + 
Vsb\ a d a + H.c) denotes the tunneling part of the Hamil- 
tonian. The current from the normal lead flowing into 
the QD is calculated by the standard Keldysh nonequi- 
librium Green's function theory, as (H = 1)113: 

/ = -4e/m J ^T N J/^G'fc) + ^(e)} (1) 

where T N/S (e) = 2irJ2 k \Vn/s\ 2 $(£ - ^N/S,k); /jv(e) is 
the Fermi distribution of the normal lead. The subscript 
"11" means taking the "11" element of the 2x2 matrix. 
The entire analysis therefore falls on to derivations of 
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retarded and Keldysh Green's functions G r and [ G < for 
the QD in the well known Nambu representationO. 

We have solved G r (e) using the equation of motion 
method. Although this method is quantitatively less 
accurate in predicting intensity of Kondo effect, it has 
been proven to. provide correct qualitative physics at low 
temperatureal3, therefore it is sufficient for the purpose 
of this work. Then we obtain the matrix form of G r (e) 
as,E_3 



e - £d - 



UA.B^ + ^A^,) 
^s°;2i + UA 2 B(H^ 21 + S s ° 21 AiSfJ) 



.(o) 

J S;12 



where So is the exact self-energy for noninteracting sys- 
tem. Then, from £< - £> = S r - S a , it is easy to deter- 
mine X = Ti r — Ti a . This ansatz has several advantages: 
(i) it is exact both in equilibrium (the bias Vb = 0) and in 
noninteracting limit (U — 0); (ii) when. A = 0, it is con- 
sistent with the original ansatz of NgU; (iii) the current 
conservation is automatically satisfied; (iv) it guarantees 
the exact relation of matrix elements Gf 2 (e) = — G^ 1 *(e), 
which in turn guarantees < dyd^ >=< djdt >*; (v) 
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where Ai(e)(i = 1, 2) and B(e) are defined as Aj(e) = [e=F 
e d T«7-2Sf-S«]- 1 ; B(e) = [l-^A^A^; 
and the self-energies S = Sjv + £s are: 
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and E£ 6 = S£ a) + S. 
expressions 15 = eti,e 
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where index i = 1,2. In these 
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Pi =< b- ki b k] >, 



and p 2 = p*. The self-energy Ejv for the coupling to 
the normal lead can easy be obtained from £g by setting 
A = and substituting (Vjs, esk, hv) by (Vn, ejvfc, a k(J ). 
The quantity n a in Eq.(|2j) is the intradot electron occu- 
pation number of state ct; < d^d^ > and < djdj > are 
the pair correlation in the QD due to the well known 
proximity effect. |— These quantities must be calculated 
self-consistentlyD't3. We emphasis that Eq.(|^) is suit- 
able for arbitrary Coulomb interaction strength U and 
superconducting gap A. In the zero-gap limit so that 
the superconducting lead becomes a normal lead, Eq.(||) 
reduces to that of the N-QD-N systeixM 

Next, we solve the Keldysh Green's function G < (e) 
which, for interacting systems, can not be obtained from 
the equation of motion without introducing additional 
assumptions. We follow the most commonly used ansatz 
for interacting-Lesser (greater) self-energy S a (a =< and 
>), due to NgJiH, but we generalize this ansatz to meso- 
scopic hybrid systems in the following way, 

S« = ~ [S Q (S5 - EgJ^X + X(SS - SoT 'Sol ( 3 ) 



Gji(e) and G 22 (e) automatically becomes purely imag- 

proposed to 
At least for 
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inary. In this regard, we note that Ref. 
use only the first term of Eq.(||) as S < 
our purpose this choice is not correct because our matri- 
ces in ([}]) are non-diagonal therefore cannot be permuted, 
and more importantly, it generates result violating points 
(iv) and (v). Finally, using (||), G < is obtained from the 
Keldysh equation G< = G r S<G a where G a = (G r ) t - 
With Green's functions G r and G < , from Eq. (jl|) the cur- 
rent / is calculated immediately. 

To gain physical insights to the analytical result, in 
the rest of the paper we discuss them numerically. In 
our numerical calculations, we assume square bands of 
width 2W so that T N/S (e) = T N/S 9{W - |e|), with 
W = 1000 > max(k B T,V b ,r,A). We emphatically in- 
vestigate the case when A > T$ and with asymmetri- 
cal barriers, Ts > Tn- In this case, an electron with 
energy |e| < A in the QD undergoes multiple Andreev 
reflections before it decays to the normal lead, and it 
cannot decay into the superconductor due to the gap. 
These multiple reflections give rise to Andreev bound 
states in the QD which are indicated by the peaks in the 
LDOSo. Fig. (la) shows LDOS at a high temperature 
T = .5. If U = 0, two Andreev bound states emerge 
at ±^e 2 +T 2 /A (when A » T s > IV). When U be- 
comes finite, due to a competition of intradot Coulomb 
interaction and Andreev tunneling, each U — Andreev 
bound state is split into two sub-states an energy U apart. 
As a result, four Andreev bound state peaks emerge in 
the LDOS each with half- width set by ■ Note that 
for large interactions U — > oo, Andreev reflections are 
Coulomb blockaded therefore only two ordinary reso- 
nance peaks at Cd and ed + U can be detected in LDOS. 

Next, we drop the temperature to T = 0.005 so that 
Kondo effect can be investigated, the data shown in 
Fig. (lb). In LDOS, the four broad peaks correspond to 
the four Andreev bound states just discussed. On top of 
them, four additional narrow peaks emerge. These nar- 
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row peaks only exist at low temperature and they have 
the typical character of Kondo effect: their peak heights 
greatly depend on T and they heighten with reducing T. 
The positions of these Kondo peaks are at e = /it at, -^jv, 
2ed + U — (Uat, and —2ed — U + [In- Note that two of 
them are dependent on a which is qualitatively differ- 
ent from the conventional Kondo effect. Moreover, the 
Kondo peaks at — (j,n and —2ed — U + hn arc outside the 
Kondo region of a normal system. Thus, instead of one 
Kondo peak, we now have three "excess Kondo peaks" . 
At bias Vb = 0.7, the two Kondo peaks at 2ed + U — fiN 
and — 2e<j — U + overlap each other, therefore only 
three peaks are observed (see Fig. (lb)). As an extremely 
strict confirmation of our analytical derivations and nu- 
merical calculations, for all cases we have checked that 
J de[LDOS(e)] — 2 which indicates that there are two 
states in the QD. 

Where do the excess Kondo peaks originate? They 
certainly cannot be due to first-order tunneling processes 
because Andreev bound states do not align with [In- In 
addition, although direct Andreev tunneling does occur, 
their effect is to evolve the intradot level td into the An- 
dreev bound states, i.e. to give the four broad peaks in 
the LDOS as discussed in the above. They cannot give 
rise to the Kondo effect simply because this process can 
not flip local spin. 

Our investigation suggests that the excess Kondo peaks 
originate from an interesting co-tunneling process not 
discovered before, which is indicated by Fig. (2). First, 
the conventional Kondo peak at e = /Ujv originates from 
the normal co-tunneling event between the QD and the 
normal lead depicted in Fig. (2a). Note that this pro- 
cess does not induce any net current. Second, the ex- 
cess Kondo peak at e = — /ijv is from a Andreev-normal 
co-tunneling process described in Fig.(2b-d). To start, 
an electron with spin-up, for example, occupies the QD 
(Fig. 2b). Then a down-spin electron in the normal lead 
with energy ~ hn can tunnel into QD and reach the 
QD-S interface causing an Andreev reflection by which 
a hole is reflected back, Fig. (2c). If Yn > Ts, this hole 
can easily tunnel into the normal lead and neutralizes an 
up-spin electron. This is a real Andreev process in which 
two electrons with opposite spin in the normal lead are 
annihilated concomitant to the creation of a Cooper pair 
in the superconductor, leaving the occupancy and the 
spin of the QD unchanged. However, when Ts > Tn 
which is our concern, it is more difficult for this hole to 
tunnel into the normal lead, the probability is increased 
for this hole to combine with the original up-spin elec- 
tron of the QD (Fig. 2c). Note that this is a virtual An- 
dreev tunneling process: after it the system is in a high- 
energy virtual state which can only exist for a timescale 
<~ U/\fj,N + €d\- Closely following this virtual Andreev 
process, another virtual process involving a normal tun- 
neling event can occur where a down-spin electron with 
energy ~ — fiN tunnels into the QD. As consequences of 
this virtual Andreev and virtual normal co-tunneling pro- 
cess, two electrons with same spin in the normal lead are 



annihilated, a Cooper pair is created in the superconduc- 
tor, the spin of the QD is flipped and a net current flows 
through the QD. At low temperatures, a coherent super- 
position of all possible co-tunncling processes of this kind 
produces the Kondo effect. The entire system, QD plus 
leads and including the Cooper pair, forms a spin sin- 
glet. The local spin is effectively screened and two very 
narrow Kondo peaks emerge at e = hn and — /it at in the 
LDOS. The new physics here is that Andreev tunneling 
can directly participate co-tunneling processes, in other 
words the Cooper pair can directly participate formation 
of the spin singlet. 

Third, the excess Kondo peak at 2ed + U — ^n is from 
the Andreev assisted co-tunneling process. An Andreev 
event involves a two-particle tunneling process, hence the 
electron number in QD can fluctuate by two. In partic- 
ular, when Ts > Tn, Andreev tunneling occurs more 
frequently than normal tunneling, therefore the proba- 
bility for QD to have zero- or two-occupancy is greatly 
enhanced. This induces the following co-tunneling pro- 
cess. When QD is empty, an electron with energy ~ ^n 
in the normal lead can tunnel into the QD state e<j, fol- 
lowing closely by another electron having opposite spin 
and energy <~ 2ej + U — ^n tunneling on to the state 
Ed + U. If the QD begins with a two-occupancy, the in- 
verse process can occur. These co-tunneling processes 
induce the Kondo peak at 2e<j + U — hn- Notice that the 
electron fluctuation is two in these co-tunneling process, 
so it only occurs under the assistance of Andreev tunnel- 
ing. Finally, the Kondo peak at —2ed~U + i^n originates 
from a combination of the processes discussed in this and 
the last paragraph. 

The above picture clearly indicates the physical process 
behind the excess Kondo peaks. Is the tunneling current 
enhanced by these excess Kondo resonances? The inset of 
Fig.(lc) plots current versus Cd at bias V& = 0.5. A finite 
bias makes the conventional Kondo peak at /Ltjv not to 
align with jis, so that this peak alone contributes weakly 
to current. At low temperatures (solid curve), our results 
suggest that the current is indeed enhanced in the region 
—U to by the excess Kondo effect because the Andreev- 
normal co-tunncling does generate net current flow. Un- 
der what conditions can we observe the excess Kondo 
resonances? (i). When > Ts, no excess Kondo peaks 
can be obtained and only the conventional Kondo peak at 
Hn emerges (inset of Fig. lb). This is because the normal 
tunneling at the N-QD interface is frequent which di- 
minishes the virtual Andreev tunneling and the two elec- 
tron fluctuation processes necessary for the excess Kondo 
peaks, as well suppresses the formation of Andreev bound 
states, (ii). When U — > oo, all excess Kondo peaks dis- 
appear, because at large U the timescale (lifetime) of the 
virtual two-occupancy state, ~ h/\U + €d — Hn\, goes to 
zero thereby greatly suppressing the Andreev tunneling 
(including the virtual process). This is the qualitative 
difference between the excess Kondo peaks discovered in 
this work and the Kondo peaks at ±A when A ~ T: 
the latter survive the U — > oo limit because their physi- 
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cal origin is from co-tunneling between the QD and the 
superconducting quasi-particle spectrum and does not in- 
volve direct participation of Andreev tunneling process, 
(iii). When /in < £d, all Kondo peaks disappear and 
four broad peaks due to Andreev bound states still exist 
as shown in Fig.(lc). Although — lin is still between a 
and €d + U, the excess Kondo peak at — hn neverthe- 
less disappears because when /Ujv < e<2, the Fermi surface 
is below ed which forbids all co-tunneling process from 
happening. 

In summary, for a N-QD-S device satisfying the condi- 
tion A > T$ > Tn, excess Kondo peaks are predicted in 
the local density of states which enhances the tunneling 
current. These excess Kondo peaks are due to an inter- 
esting co-tunneling process involving an Andreev tunnel- 
ing from the QD-S interface and a normal tunneling from 
the N-QD interface. This new co-tunneling phenomenon 
clearly demonstrates the rich many-body physics of the 
hybrid system in which the formation of a spin singlet 
not only can be due to the many free electrons, but also 
can be due to a combined effect of free electrons and the 
many Cooper pairs in the superconducting lead. 
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FIG. 1. (a) LDOS vs e at different U for a high tem- 
perature T = 0.5, with parameters VJ, = 0, Tjv = 0.1. (b) 
and (c) are LDOS at different bias V& for a low temperature 
T = 0.005, and parameters U = 0.7, T N = 0.15. Other pa- 
rameters are e d = 0, Tg = 1, A = 10. The inset of (b) shows 
LDOS vs e for symmetric barriers with Tn — Ts = 0.15. The 
inset in (c) shows current vs at T — 0.5 (dotted curve) and 
T — 0.005 (solid curve). Other parameters of insets are same 
as those of the solid curve in (b). 

FIG. 2. Schematic plots showing the co-tunneling process. 

(a) Conventional normal co-tunneling in the N-QD interface. 

(b) , (c), and (d) are the Andreev- normal co-tunneling, (b) 
the initial state; (c) the virtual tunneling event of the An- 
dreev-normal co-tunneling; (d) the final state. The dotted 
arrow in the QD of (c) represents a hole, other arrows repre- 
sent electron. 
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